Spin-orbit induced chirality of Andreev states in Josephson junctions 
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We study Josephson junctions (JJs) in which the region between the two superconductors is a multichannel 
system with Rashba spin-orbit coupling (SOC) where a barrier or a quantum point contact (QPC) is present. 
These systems might present unconventional Josephson effects such as Josephson currents for zero phase differ- 
ence or critical currents that depend on the current direction. Here, we discuss how the spin polarizing properties 
of the system in the normal state affect the spin characteristic of the Andreev bound states inside the junction. 
This results in a strong correlation between the spin of the Andreev states and the direction in which they trans- 
port Cooper pairs. While the current-phase relation for the JJ at zero magnetic field is qualitatively unchanged 
by SOC, in the presence of a weak magnetic field a strongly anisotropic behavior and the mentioned anomalous 
Josephson effects follow. We show that the situation is not restricted to barriers based on constrictions such 
as QPCs and should generically arise if in the normal system the direction of the carrier's spin is linked to its 
direction of motion. 

PACS numbers: 03.75.Lm,72.25.Dc,71.70.Ej 



I. INTRODUCTION 

The study of spin-orbit coupling (SOC) and related effects 
is one of the most active fields in mesoscopic physics. As 
the fabrication techniques improve, the link between the or- 
bital and the spin degrees of freedom can be engineered al- 
lowing for the exploration of quantum phenomena at a deeper 
level.' ^ Moreover, the importance of SOC goes far beyond 
mesoscopic device design as, for example, it has led to the 
identification of the Quantum Spin Hall phase in topological 
insulators,^ a new phase of matter. ''^ More recently, hybrid 
systems involving SOC and superconductivity have become 
one of the most investigated platforms.'^ When in combina- 
tion with an appropriate magnetic field (generated externally 
or induced by proximity with a magnet in the sample) the sys- 
tem is expected to develop topological superconductivity.""'"' 
This opens the possibility to synthesize Majorana fermions in 
semiconducting quantum wires, cold atoms or devices based 
on topological insulators. 

In this work we theoretically investigate the DC Josephson 
effect in a hybrid system having superconductivity and spin- 
orbit coupling; more precisely, two s-wave superconducting 
leads (S) coupled by a two-dimensional normal region (N) that 
has Rashba spin-orbit coupling. For superconductor-normal- 
superconductor (S-N-S) junctions, it has been recognized that 
the interplay between external magnetic fields — or, alterna- 
tively, intrinsic magnetism — and the spin-orbit coupling can 
change the relation between the supercurrent and the phase 
difference between the two superconductors.Q^JI^ 

The key ingredient in the S-N-S junction studied here is that 
the normal region contains a gate-voltage controlled barrier, 
see, for instance, the Josephson junction measured in Ref 121] 
In Ref. ^72} a weak external magnetic field, applied in a spe- 
cific in-plane direction, was predicted to generate a control- 
lable phase shift in the cuiTent-phase relation (CPR), i.e., the 
system develops a supercurrent for zero phase difference be- 



tween the superconductors. As the gate voltage is reduced, 
and the barrier allows the transport of a few (more than one) 
transmission channels, the Josephson current becomes non- 
symmetrical as a function of the phase. The maximum dissi- 
pationless current, the critical current, depends on the direc- 
tion of the current and therefore there exists a current window 
on which the supercuiTent is rectified. 

In the long Josephson junction limit with transparent S-N 
interfaces, a normal region with a one-dimensional (ID) spin- 
orbit coupling combined with a magnetic field, either parallel 
or perpendicular to the SOC axis, can only produce transitions 
between Q-n junctions as it is the case for Superconductor- 
Ferromagnet-Superconductor (S-F-S) junctions!?-^'-' In this 
limit this can be understood using Kulik theory, but the 
same result follows by virtue of symmetry relations of the 
Bogoliubov-de Gennes (BdG) Hamiltonian,^^ irrespective 
of the geometry of the junction for any generic ID SOC cou- 
pling. On the other hand, the magnetic field triggers shifts in 
the CPR by incorpora ting t he mixing with other subbands in 
an effective ID model.l^^Ell 

In the junction studied here the anomalous effects arise be- 
cause the barrier in the normal region behaves as an uncon- 
ventional spin polarizer. As it was reported by Eto et al. in 
Ref 28 , a quantum point contact (or baiTier) in a multichannel 
normal material with Rashba SOC spin-polarizes the current 
without the need for magnetic fields or ferromagnetS-I^SJEIlThe 
larger the spin-orbit strength and the smoother the barrier, the 
more the polarization grows. Due to time-reversal symme- 
try, the favored spin projection reverses if the direction of the 
transport is inverted. 

At zero magnetic field, we show how these gate-voltage 
controllable spin polarizing features in the normal device in- 
duce a correlation between the velocity'''' and the spin of An- 
dreev states. As a result, the application of a magnetic field, 
along the direction in which the SOC polarizes the current, 
triggers the mentioned anomalous Josephson effects. 
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We present a simple WKB picture for the transport through 
the barrier that successfully explains the main features of the 
exact numerical results for the Andreev states including their 
spin texture. We also show that the anomalous Josephson ef- 
fects are not restricted to QPCs, they also appear in wider bar- 
riers (stripe-shaped samples) for which the transverse modes 
are so close in energy that the plateaus of quantized conduc- 
tance can not be resolved when changing the gate voltage. 
For both types of barriers, N-S junctions with reduced trans- 
parency make the normal device behave as an Fabry-Perot 
interferometer;^"* the total transmission and the polarization 
are maximized at particular values of the gate voltage that 
controls the barrier. This resonant behavior translates into 
stronger Josephson anomalies that can be tuned by changing 
the gate voltage in a small amount. 

The paper is organized as follows, in Sec. |ll]we introduce 
the model and we outline the polarizing mechanism of the nor- 



mal device, in Sec. Ill we focus on how the spin properties of 
the Andreev states result in the anomalous Josephson effects 



mentioned above, and conclusions are presented in Sec IV 



II. 



THE MODEL AND THE PROPERTIES OF THE 
NORMAL DEVICE 

A. Model for the junction 



Our model describes a Josephson junction with a QPC or 
barrier in a normal region that has spin-orbit coupling. The 
device, illustrated at the top of Figure [T] includes two su- 
perconducting leads — Left (L) and Right {R) — and a central 
normal part (N) made out of a 2DEG in which the barrier is 
created with the help of electrostatic gates. For convenience 
we assume that the three regions are restricted to a stripe ge- 
ometry, i.e. the potential becomes infinite outside the region 
< y < Wy. The Hamiltonian describing the leads is 



/"-{ 



(1) 



Here 'y-L,R and the integral is done in the x < xi semi-space 
or the X > xk semi-space for the left and right lead, respec- 
tively. In the above expression, is the electron mass in the 
superconductor, yu is the chemical potential and "^Irif) creates 
an electron with spin cr at position r - {x,y). We assume that 
within the leads the superconducting order parameter is con- 
stant and set A^ir) - Ay and AR(r) = Aqc"^. In the following 
we take Aq - 1 .5meV as in Nb. 

The normal part of the junction, xl < x < xr, is described 
by the Hamiltonian of a 2DEG with Rashba SOC and the po- 
tential y(r) that defines the barrier, 

Hn^ \ drY^<lfl{r)'HN'iAr), (2) 

^ a-,iT' 



where m* is the effective mass of the 2DEG carriers, a is the 
strength of the Rashba coupling and cr, is the /-component of 
the spin operator. The potential V{r) includes a constant term 
that shifts the bottom of the conduction band of the 2DEG 
with respect to the one in the superconductor. 

As in Ref 22 we consider thin superconducting films paral- 
lel to the (x,3')-plane. In this configuration in-plane magnetic 
fields lead only to Zeeman interaction. Therefore, the contri- 
bution to the Hamiltonian due to an external in-plane magnetic 
field is 

Hz^ r^/r2*>)^^(B,(r, + B,,(r,,)*^,(r), (3) 

o-.a-' ^ 

where B, are the components of the magnetic field and the 
gyromagnetic factor 



gN, for xi< x< xr 
gs , otherwise 



(4) 



is different for the leads and the normal region. In the regime 
of parameters studied in this paper the superconducting con- 
tacts are almost unaffected by the magnetic field. Namely, 
we take the Zeeman energy Ez on the normal region, always 
smaller than 20% of the superconducting gap Aq in the S con- 
tacts. As the g-factor in the In-based 2DEG is at least 5 times 
larger that the g-factor in the superconductor, the Zeeman en- 
ergy in the superconductor is at most Ao/25. 

In the present work we solve the full Hamiltonian, H-Hz+ 
Hn+Hl+Hr, numerically using a finite difference approach.!^ 
This is done by writing a tight-binding version of the Hamilto- 
nian and using a recursive Green function method. The de- 
tails of the approach are presented in Appendix [A] For Aq 
one must work in a Nambu space representation since H,oi is 
indeed a BdG Hamiltonian.^^ 

In the lattice model the tunneling between the normal and 
superconducting materials is described by hopping Hamilto- 
nians, H^l and Hj^^r. The sites at the interfaces are connected 
by a spin-conserving matrix element, f/,, that arises due to the 
kinetic term in the original Hamiltonian (oc p^). As shown 
in Eq.( |A7| i, there are also spin-flipping tunneling elements be- 
tween those sites. These arise due to the symmetrization made 
to assure Hermiticity of the SOC termP^l 



a(x) 



-PxO-y ^ - r^r iaix)pj, + pM^)) O-y (5) 

2n 



Clearly, such symmetrization is relevant only where d^a(x) + 
0: here the interfaces between the normal region and the su- 
perconductors. 

As assumed in the classic Blonder- Tinkham-Klapwijk 
(BTK) approach for studying Andreev reflection,E2 non-ideal 
transparencies can be introduced through a sharp potential at 
the N-S interfaces in the continuous model. In the lattice 
model, the transparency of the barrier is controlled by the pa- 
rameter A such that 



(6) 
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where ts and tf4 are the hopping parameters for the supercon- 
ducting and normal regions, respectively. Notice that for real- 
istic junctions and f^v are different since they are chosen to 
simulate the Fermi velocity and electron mass mismatches be- 
tween the 2DEG and the superconductor A transparent junc- 
tion corresponds to A = 1 since in such a case tb is the average 
of the hopping in the two regions minimizing the scattering at 
the S-N interfaces; choosing A < 1 allows for the simulation 
of poorer transparencies. To achieve the full transparent junc- 
tion case (Z = case of BTK) it is not enough to set A = 1 
but we also need to assume that - 1]^. Besides the interface 
parameter A, the junction is characterized geometrically by its 
total length, Lf^-XR- xi, by the properties of the barrier (see 
Appendix [a]) and by the number of transverse channels in the 
normal region. 

B. Normal device: spin polarization due to a barrier potential 

Understanding normal transport through smooth barriers, in 
presence of spin-orbit coupling, is a prerequisite to the inves- 
tigation of the JJs we treat here. Therefore, in this section, we 
focus on the case in which the leads are not superconducting 
but instead they are: (/) a continuation of the normal region, or 
(//) a metallic material M, different from the spin-orbit cou- 
pled 2DEG in the central region. For the case (/) Eto et al?^ 
have shown that a 2DEG with a QPC is able to spin polar- 
ize the current without the help of an external Zeeman field. 
They assume a Rashba SOC, and study the ballistic transport 
through QPCs described by a quite general saddle-like poten- 
tial. 



1. Origin of the current polarization: a qualitative picture 

We focus on case (/) for which the interfaces 2DEG-leads 
do not introduce backscattering. The main underlying mech- 
anism can be exemplified using two conducting channels n - 
1, 2 in a Wentzel-Kramers-Brillouin (WKB) picture.!^ At ev- 
ery position Xa one obtains the dispersion relation that would 
follow if the system were infinitely long with the y-dependent 
potential, V^Sy) = y(xa,y)- We first exclude from the Hamil- 
tonian of the normal region in Eq.Q the part of the Rashba 
SOC proportional to the perpendicular momentum, 'Hf - 

fPyO-_,. For "K" = ^ - fpxO-y + Vjc^^iy), solutions at x„ are 
uniquely identified by k^, by the spin projection along the y- 
axis, cTy - +, and by the transverse mode number n. For each 
transmission channel, (n, o",,), the kj^ solution at position x^ is 
\'^x„n}\kt}\o-y), with eigenenergies, 

El.a/^a) = ^ (k, - k.a,)' -E, + E^(xa) , (7) 

where ka = = - ^0- and E^(xa) are the trans- 

verse mode energies, namely, the eigenenergies corresponding 

to the eigenfunctions, (ylO j^,,,), of the Hamiltonian ^+Vx„(y)- 
Importantly, the E^^^^ ^ ^(xa) parabolas for modes (1, -) and 
(2, +) (as well as for (1, H-) and (2, -)) cross each other; see 



gray lines in Fig. [TJa). Away from the center of the bar- 
rier and for strong enough Rashba interaction these cross- 
ings are below the Fermi level, Ep. The transverse SOC, 
'H'^ oc pyCTx, produces avoided crossings because it coher- 
ently mixes transverse modes with different parity, n mod 2, 
(since p,, = and opposite cr,.. Figure 1(a) sketches the 

transport through the QPC in the adiabatic picture. For the 
situation shown there, the value of the potential at the sad- 
dle point, Vg, is tuned to allow the passage of a single trans- 
verse mode: the conductance is 2e^/h. As electrons at Ep 
pass through the QPC from left to right their k^ momentum 
change. Electrons in the lowest-lying split channel {!,+) pass 
through unaffected. On the other hand, electrons in the higher 
split channel (2, +) are swept twice through the anticrossing: 
first changing to (1, -) and then to (2, +). Incoming electrons 
occupying the down modes (1, -) and (2, -) (see contact with 
the Left lead at position xl) do not arrive at Xc (center of the 
barrier) but instead they are backscattered to xl into the left- 
traveling modes (2, +) and (1, -f), respectively. Thus, the cur- 
rent arriving at xr (the position of the contact with the Right 
lead) is up-spin polarized along the y-direction. 

In the two-modes sketch of Fig. [T]we have assumed that 
the avoided crossings are fully effective. In a more general 
case the two quantum channels associated to the first trans- 
verse mode can be divided in: (i) channel (1,-1-) which is fully 
transmitted to xr without state conversion, and (ii) combina- 
tion of (2, +) and (1, -) that arrives to xr after mixing in region 
of the barrier. The larger the strength of the Rashba coupling 
and/or the smoother the barrier potential, the more effective 
are the avoided crossings and the better becomes the polarizer: 
the contribution to the current of the (2,+) quantum channel 
is increased. As the polarization is defined by the transport 
properties of electrons at the Fermi level, what matters most 
is the gradient of the barrier potential at the positions x^'^ in 
which the crossings between the relevant dispersions fall at 
Ep, i.e.. 



^F = 4,^„)(^^=) = £[U^,)(x-). (8) 



Two positions, x'j^ < Xc and x'^ > Xc, fulfill this condition. 
If the barrier is symmetric with respect to x,. both avoided 
crossings share the same effectiveness. In Ref. |28]the prob- 
ability of state conversion is estimated — using the Landau- 
Zener formula for the case that Vxiy) is a hard-wall constric- 
tion potential of variable width W{x) — as Plr-^— e^^"'^ with 
^*A:„|H^(4^)/[f (4y]|. 

Notice that for barriers with coexistence of Rashba and 
Dresselhaus SOCs, nonzero polarization arises in exactly the 
same way along the direction of the spin operator that mul- 
tiplies pg, with = (cos 6, sin ff), the direction parallel to the 
current that flows through the barrier (here 6-xand the corre- 
sponding spin operator for the Rashba SOC is o-y). This opens 
the possibility to electrically change the spin-polarization axis 
of an Eto et al. polarizer as the Rashba strength is tuned with 
gate-voltage'*' and the Dresselhaus interaction remains con- 
stant. 
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FIG. 1. (Color online) (a) UPPER PANEL: scheme of the junction including the leads and a qualitative color plot of the barrier potential, 
V(r); Vg is the value at the saddle-point. LOWER PANEL: qualitative two-mode picture of how the current through the barrier becomes spin 
polarized (left to right case). In the dispersion relations E{k;,) at different positions we indicate the relevant states using dots at the Fermi energy, 
iSp. Electrons from the Left lead (at x = xl) are injected in four states, (n, CTy), with « = 1, 2 and cr,, = ±. The barrier allows only one transverse 
mode to be transmitted (first conductance plateau, G = 2e-/h). Electrons injected atx = xi in states (I, -I-) and (2, -I-) are transmitted to the right. 
They arrive at x = xr in states (1, +) and (2, -I-) and therefore the current is up polarized along y-direction. Notice that — in the passage through 
the barrier — an electron injected in state (1,-1-) does not see any avoided crossing (as in a zero-SOC barrier) whereas an electron injected in 
state (2, +) suffers two SOC-induced avoided crossings: first it is mixed with (I, -) and then it returns to (2, +). (b) Summary of the left to 
right and right to left transport properties of the barrier. The spin polarization from right to left is down as electrons injected (at x = xr from 
the Right lead) in (1, -) and (2, -) are transmitted whereas electrons injected in (1, +) and (2,-f) are reflected to (2, -) and (1, -), respectively. 



2. Time-reversal symmetry and transport from right to left 

In the above discussion we have focussed on the left (L) to 
right (R) transport in the symmetric case. We now relax the 
condition of spatial symmetry of the barrier profile allowing it 
to be different for x < and for x> Xc. Therefore, we assume 
that the probability to follow the avoided crossing is p/, at X^^ 
and pj{ at id^. For the two-channel model, with the barrier 
in the first plateau of conductance, the nonzero transmission 
coefficients (from left to right) at the Fermi energy are 

(9a) 
(9b) 



1£l 

h ' 



Pl-^r = PR and 



of Eq.(j9| in Eq.([TO]l obtaining: Gl_ 

Importantly, due to time-reversal symmetry the resolved 
transmission coefficients at the Fermi energy fulfill the fol- 
lowing relation 



(L,«t,iTz,)->(«,«s,o-«) 



(«,nR,o-fi)->(L,«i,o-i) 



(11) 



7"(L,1,+)^(R,1,+) - 1 



{L,\-y 



AR,\-) = (1 - ft)(l - Pr) 



(9c) 
(9d) 
(9e) 



r(L,l,-)^(R,2,+) = (1 - Pl)PR ■ 

T(L.2,+)^(R,i-) - PLii " Pr) ■ 

T(L,2,+)->(R,2,+) - PlPr ■ 

The conductance, G, the polarization of the transmitted elec- 
trons, P, and the relative spin polarization of the incoming 
electrons that contribute to the current, D, can be expressed in 
terms of the transmission coefficients as, 

,2 



e 



' (L,ni,(ry^L)^(R,nK,CTy,,) 



2 



hG 



nL,a-yL)^(R,n]t,CTy_f:) 



.nL,o-.vx)->(«.««,o",,«) 



(10a) 



(10b) 



(10c) 



where the subscript L ^ R refers to the left to right trans- 
port. For the polarizing barrier, we introduce the coefficients 



Therefore the resolved transmission in the opposite direction 
of electrons with opposite spins (at both ends) is identical. 
This follows from the self-duality of the scattering matrix in a 
spin- 1/2 time-reversal symmetric system.'^ Here, this expres- 
sion is of great importance for the Josephson effect because 
at the S-N interfaces an electron coming from the normal re- 
gion can be Andreev-reflected back to the normal as the lack 
of an electron with opposite spin (a hole) in the same trans- 
verse mode. Since this hole propagates as the lacking electron, 
Eq.([TT| establishes the main difference with a conventional 
magnetism-based polarizer In the latter systems, irrespective 
of the direction of the current, the largest transmission is for 
electrons with spins parallel to the majority-spin. 

For the polarizing barrier treated here, by combining the 
time reversal relation with Eqs.(j9]l and ([TOjl, we get that the 
normal transport from right to left has Gr-^l - Pr-^l - 
-pi and Dr^l - -pR. Note that for the case of the smooth 
symmetrical barrier with strong Rashba of Fig. [T]/7l - Pr~^ 
and thus Pr^l - ~Pl-^r ~ 1- The inversion of the sign of 
the polarization is summarized in Fig. [ijb). It can also be 
understood from the WKB picture in Fig. [TJa) by noting that 
the electrons arriving to L are due to incoming electrons from 
R entering the sample: (i) in channel (1,-) passing through 
the barrier without mixing or (ii) in channel (2, -) arriving to 
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L in channel (2, -) after mixing in the region of the barrier 
with channel (1, +). 

We have introduced the two most important differences be- 
tween the Eto et al. polarizer and a magnetism-based spin 
polarizer. First, here the device polarizes even though the con- 
ductance is equivalent to two quantum channels, this contra- 
dicts the naive expectation that G-2e^/h implies spin degen- 
eracy and that therefore the polarization should be zero. The 
second distinctive feature is that in symmetrical barriers the 
spin polarization of the transmitted current reverses when the 
current flow is reversed. More generally, irrespective of the 
spatial symmetry of the barrier, we have that Pl^r = -Dr^i 
and Pr^l - -Dl^r which means that the SOC induced po- 
larization does not conspire against the formation of Andreev 
states. The preferred spin direction for the transmitted elec- 
trons arriving to a given N-S interface is exactly opposite to 
the preferred spin direction for the incoming electrons (or lack 
of them) that can be transmitted to the other side of the barrier 



In Sec. Ill we study which properties of Andreev states are af- 
fected when these SOC-based polarizer barriers are present 
inside a S-N-S Josephson junction and how these properties 
manifest when a weak Zeeman energy is present. 



3. Numerical results in more realistic models 

Using the QPC potential of Eq.( |Al| l and the tight-binding 
model we compute the conductance and the polarization in 



devices with many channels. In Fig.|2]results are presented as 
a function of Vg for case (/) and (//) described above. We have 
chosen the opposite limit to the one shown in Fig. [T[ the SOC 
strength is weak, the avoided crossings are not fully effective 
and thus the polarization is not strong. Case (/) is shown in 
Fig. |2|a), as the leads are 2DEG-based the scattering back to 
the barrier is minimized and the plateaux of conductance are 
well defined. The SOC-induced current polarization is always 
smaller than 5%, the characteristic value in the first plateau is 
P('/)' = 0.04. 



Figure[2|b) shows results for the case (//) in which the leads 
are metals and therefore the electrons at the 2DEG-M inter- 
faces are strongly scattered back to the 2DEG region. Reso- 
nances in the conductance are observed for special values of 
Vg. Correlatively, the spin polarization of the current shows 
maxima too. Notice that the average spacing between these 
maxima is larger than that of the conductance maxima. The 
main result is that constructive interferences, due to multiple 
reflections of channels between the interfaces and the QPC, 
strongly boost the current spin polarization with respect to 
case (/). We have checked (not shown) that the relative po- 
larization enhancement, P^^jj^ IP\i\, grows the weaker is the 
SOC strength. It is interesting to analyze this amplification 
mechanism for the same phenomenon is correlated with the 
anomalous Josephson current amplitude in the S-QPC-S case. 



(a) case (I) 
2DEG-QPC-2DEG 




V (meV) 

g 



V (meV) 



FIG. 2. (Color online) Conductance G (gray line) and polarization P 
(black line) as a function of Vg for the quantum point contact, QPC2, 
described by Eq.([ATJ taking {z\Wh;Lt) = (90; 100; 250)nm. The 
2DEG parameters are for an InAs-based 2DEG, (n = lO'^cm"^ and 
m* = 0.045mj with m^ the electron mass) Ep = 53.3meV, a = 5meVnm 
and L/^ = xjf-xi^ = 1 .2yum. In the simulation we take IV,, = 1 80nm and 
flo = 3nm. (a) Case (/), 2DEG-QPC-2DEG: the leads are the same 
2DEG and therefore the scattering at the interfaces is absent. As the 
Rashba strength is weak the polarization is small: less than 4% in the 
first plateau, (b) Case (//), M-QPC-M : The leads are a metallic ma- 
terial, M, different from the 2DEG leading to a strong backscattering 
at the M-2DEG junctions; taking metallic Nb as a reference we use 
ra*/m, = 0.045 and the Fermi velocities ratio, Vp'-'^^/Vp'' » 2.25. Due 
to Fabry-Perot interference-like effects both the conductance and the 
pofarization show resonant behavior as a function of V^; remark- 
ably, there are values of Vg in which the polarization can be much 
larger than the sub 4% values observed when the device is in between 
2DEG leads. 



A simplified scattering theory is sufficient for a qualita- 
tive understanding.^'* It includes the two lowest channels, the 
M-2DEG interfaces with a reflection coefficient ro, the SOC 
strength through the probability p of state conversion at the 
avoided crossing between channels ( 1 , cr,,) and (2, iXy), and the 
phases accumulated through the QPC for each of the quan- 
tum channels. These phases can be calculated by a WKB for- 
mula, they depend both on the junction length and on the gate 
voltage. The result of this calculation is that the full trans- 
mission (thus the conductance) as well as the spin polariza- 
tion oscillate with the accumulated phases due to Fabry-Perot 
interferometry-like effects. As in the exact full numerical sim- 
ulation presented in Fig. [2] for a much larger total number of 
channels, taking a large ro value one sees: (i) the ratio be- 
tween the polarization value at the peaks and the polarization 
for ro = is larger the weaker is the Rashba strength; and 
(ii) the oscillation period for the polarization is larger than 
the one of the conductance. The explanation for this is that 
the spin polarization of the current crucially depends on a co- 
herent mixing of split channels with opposite spin directions. 
The polarization is controlled by the phase accumulated in the 
region between the output interface, at xr, and the closest mix- 
ing location, at x'^^. On the other hand, the full transmission is 
sensitive to the phase acquired through the whole length of the 
device, thus explaining the faster variation of the conductance 
with the gate voltage.^ 
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III. THE EFFECT OF SOC-BASED POLARIZERS IN 
JOSEPHSON JUNCTIONS 

Defining the superconducting phase difference, (f> = (pi- (pjf, 
the Josephson current is given byj^^EH 



J{e,<p) 



2e IdH 



J{s, (p)Me) de ., 



(12) 



where (...) indicates the expectation value at thermal equilib- 
rium , fti(x) =1/(1 +e''^**^) is the Fermi distribution and we 
have introduced the Josephson current density 



ds 



7(e,0) = 2-^(5(£-e,(<^)). 



(13) 



Thus J(s,<p) contains the information about the velocity at 
which Cooper pairs (i.e., charge 2e) are transmitted via An- 
dreev bound states, namely, the states with energy Sj which 

have nonzero Ij^- As we discuss below, the eigenvalues ej of 
the system do not need to be computed in order to get the An- 
dreev levels. It is sufficient to compute J{s, </>) and later iden- 
tify the coordinates (e, 0) at which J(s, (fi) contributes to the 
current. In the discretized model presented in Appendix [A] the 
current is obtained from Eq.( |Al \\ . The calculation involves 
expectation values that can be expressed as integrals over e of 
retarded and advanced Green function.'*"'. This method allows 
us to directly obtain J(b, <p) and from it identify the contribut- 
ing Andreev energy levels. 

Notice that in the regime of parameters we work, 
47r-periodic components are not present as the magnetic field 
is weak and it is applied parallel to the direction of the 
SOC field; this is not the regime appropriate for Majorana 
physics !^ ' I ^'^ l 

Here, due to the SOC -polarizing barriers and by virtue of a 
small magnetic field along the polarization direction, we find 
strongly asymmetric (2;7r-periodic) I{(p) shapes. We have veri- 
fied that the shapes always fulfill 



Jo ' 



I{4>)d(f> ■■ 



2e 



{{H)^=2.-{H)^=o)^0, (14) 



i.e., the area below the CPR is zero. This trivial property is 
not evident from the plots reported in the following sections. 



A. Anomalous Josephson effects enhanced by 
barriers with SOC: ideal interfaces 

In this section we study the case of ideal S-N interfaces 
in order to isolate the effect in the Josephson current due to 
the presence of SOC-induced spin polarization — in Sec IIIB 



we will focus on how the basic properties discussed here are 
modified by the resonant-like behavior that appears for more 
realistic interfaces.i^ The effective mass and Fermi velocity 
of the superconducting leads are matched to the ones of the 
2DEG and normal reflection (for |e| < Aq) is minimized. In 
the lattice model this is achieved by setting -ti^-ti, [A = l 
in Eq.(|6]l]. 




Josephson curren t 



/] 


(c) 
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s/\s I 
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15 




(^71 1 X (a^) 

FIG. 3. (Color online) Exact results for the full multichannel model 
of a Josephson junction with a SOC-based spin-polarizing barrier. 
The superconductors have Aq = 1.5meV and are chosen to minimize 
normal reflection (see text). The normal region length is Ljv = 300ao = 
900nm and the QPC parameters are those given in Fig. [2] Here, the 
gate voltage is set in the middle of the first conductance plateau and 
the SOC is strong, a = 40meVnm, so that the polarization is close to 
one. Color maps as a function of energy, e, and the superconducting 
phase difference if> for: (a) Px^is, (p), the local density of states at the 
last layer of the normal, Xn; (b) and (d) /(e, 0), the current density. 
Notice that J(e, 0), for |£| < Aq, images the Andreev states: those that 
contribute to the current (red and blue — in the online version — are 
positive and negative contributions, respectively, given by the sign 

d£ 

of the slope -^). (c) Total Josephson current, /, i.e., the integral of 

J{e,(I>). Panels (e) and (f): magnetization, Sy(x,y) = |<it,), at 0o = 
03n for Andreev states |A) and |6), (see panel (d) and the discussion 
in Fig.|4]l. The magnetizations, integrated over all the normal sample, 
are positive: = 0.45 i^. 



1. Spin properties of Andreev states in absence of magnetic field 

Figure |3] shows the results for zero magnetic field for a de- 
vice with strong SO coupling. In Figj3|a) we show a map of 
the local density of states (LDOS) at x-xr as a function of e 
and (p: px,i(s,(p) (see Appendix [A|. The non-dispersive states 
are due to modes that are reflected at the barrier and do not 
contribute to the Josephson current density shown in Fig[3jb). 
For the gate-voltage shown, the barrier allows the transmis- 
sion of a single transverse mode. Thus, at each side of the 
barrier there are standing waves that can not transport Cooper 
pairs from one superconductor to the other — except for an ex- 
ponentially small contribution due to quantum tunneling. 

On the other hand, the Andreev bound states — those con- 
tributing to the current — have an apparent linear behavior with 
(p. The sign of their contributions to J(s, 0) is correlated with 
the sign of their slopes: the latter are linked to the veloc- 
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FIG. 4. (Color online) Sketch for the WKB interpretation of the Andreev states |A) and \B) shown in Fig.|3[d), (e) and (f). Both states have 
positive velocity, they transfer Cooper pairs from the Left (L) to the Right (r) superconductor. Figure |3] shows that the state |A) transports 
Cooper pairs at a faster rate (since ^ > ^) and its Sy profile is always positive with a single maximum along the j-direction, whereas Sy for 
the Andreev state |S) shows two positive maxima except from the center of the barrier where there exists a negative contribution to Sy having a 
single maximum. This is in full agreement with the WKB picture presented here. For the state |A> (left panel) the electron occupies the (1, +) 
transverse mode all along the barrier as much as the reflected hole occupies the (1, -) transverse mode. On the other hand for |B> (right panel), 
because of the smoothness of the barrier and the avoided crossings presented in Fig.[T] the electron (hole) is on the (2, +) ((2, -)) mode away 
from the barrier and it mixes with the (1, -) ((1, +)) close to the barrier. This explains the different patterns of local magnetization along the 

electron-to-hole Andreev reflection at R, [Hi) hole traveling 
from R to L, and (iv) hole-to-electron Andreev reflection at L. 

For the subgap Andreev transport (|e| < Aq), we sketch in 
Fig|4]the two diff'erent possibilities left open by the device as 
Vg is tuned to allow a single transverse mode to transmit. The 
left panel shows different processes that lead to the state \A): 
the electron travels from L to R in the (1, H-) mode and it is 
Andreev reflected as the lack of an electron in the (1 , -) chan- 
nel; then, this hole is transmitted in the (1, -) channel from R 
to L where it is Andreev reflected as an electron back in chan- 
nel (1,-1-) and the loop closes. This is therefore an state with 
an overall spin-up polarization in the y-direction and a trans- 
verse spatial profile given by the first transverse mode (which 
has a single maximum in the transverse direction). This is in 
agreement with the exact results for the spin profile shown in 
Figige). 

Similarly, the right panel of Fig|4]focuses on the case of the 
Andreev state \B}; the electron in its travel from Lto R starts 
at mode (2, +), it mixes with the (1, -) mode in the region of 
the bamer and then arrives at R in the (2, 4-) mode. It is An- 
dreev reflected as the lack of an electron in (2, -), this hole 
is transmitted arriving at L in (2, -) but in the region of the 
barrier it mixes with mode (1,+). This is also in agreement 
with the exact quantum result shown in FigjSj^f) for state \B}, 
as the local spin density shows two positive maxima (spin up 
and second transverse mode) away from the barrier and a con- 
tribution with a single negative maximum (spin down and first 
transverse mode) appear in the region of the barrier. Further- 
more, within this picture, it is also justified that the Andreev 
state |A) transports Cooper pairs at a faster rate than \B) [see 
greater slope with (p in Fig.[3]^d)] because in the first one both 
electrons and holes are transmitted fully in the first transverse 



y-direction, and why the Andreev state |6) is slower than the |A) one. 

ity of the Andreev states [cf. Eq. ( [T3| ), see Figj3|b)]. In 
Fig|3jc) we present the total Josephson current, /, which is 
given by the integral in Eq.([T2]i. It includes the contributions 
of both the discrete (Andreev bound states) and of the con- 
tinuum (|e| > Aq) spectrum. The most impo rtant contribution 
comes from the Andreev bound states. ^^"^^ The shape of the 
CPR for zero field does not show signatures of the presence 
of the SOC and the baiTier: similar triangular-like CPR shapes 
are expected in S-N -S long ju nctions with high transparency 



at the S-N interfaces.EHEUt 



Even though the magnetic field is zero, the Andreev levels 
are not spin-degenerate [Figj3|d)]. This is due to the Rashba 
coupling in the normal region and — as it becomes clearer 
below — the effect is enhanced by the presence of the spin po- 
larizer barrier. Let us fix = O.Stt and investigate the features 
of two neighbors Andreev states, |A) and \B}. Notice that 
both Andreev levels have the same slope sign and therefore 
they transport Cooper pairs in the same direction. In a sys- 
tem without magnetism and without SOC, the two states |A) 
and \B} would be degenerate and their spin properties would 
be opposite to each other. Here the situation is different: the 
two Andreev states neither share the same transversal profile 
nor their spin-properties are opposite: in Figs.|3je) and[3jf) 
we show their spin density along the y-direction — the direc- 
tion in which the barrier polarizes the current as a result of the 
Rashba coupling. 

In order to qualitatively analyze the results it is useful to re- 



sort to the WKB picture presented in Sec II B First, we notice 
that in the states \A) and \B) the Cooper pairs are transferred 
from the Left (L) superconductor to Right (R) superconductor 
(positive slope). Thus, each Andreev bound state is a station- 
ary cyclic sequence of: (/) electron traveling from L to R, (ii) 
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mode, which is at all positions faster than the other channels, 
as it is evident from the slopes I ^| in Figjl] 

2. Symmetry considerations and the weak magnetic field limit 

We start studying the symmetry properties of the BdG 
Hamiltonian, Hbcic, in absence of ferromagnetic materials and 
magnetic fields. In this case, we have 



(15) 



where T is the time-reversal operator Notice that -(p is equiv- 
alent to 2n — (p. With the exceptions of = and (p-n, Eq.([T5|l 
does not imply the existence of degenerate solutions for Hsdc 
at a fixed 0. Instead, it allows us to construct solutions for the 
case of = (In - (p^) from the ones at = 0^. Thus, for any 
given eigenstate (labeled by the natural number j-1,2, ...) \ j} 
at (e, 0) = {sj, (px), we know that 



I;) 



\])=T-\i) 



(16) 



This means that if the energy of \i) aX. 4> - S4> + (p.x is given 
by Ej{5(p) - Sj{(px) + bj 6(p then the energy of | j) at (p = S(p + 
(In - (px) is Sj(64>) - Sj((px) - bj 5(p, with bj the slope of the 
Andreev level for |j) at = (p^- Therefore, their velocities, Vj 
and VJ, become opposite to each other (v^- - -vj) as they are 

proportional to ^ and respectively. This implies that the 
Josephson current must be zero for either 0, - Q or (p - n, as 
the two 7~-related states share the same coordinates (e, (p) so 
their individual contributions cancel out. For arbitrary values 
of (f) the following relations hold 

/(0) = -/(-0) = -/(27r-0), (17a) 
I{n + (p)^-I{n-(p). (17b) 

An obvious property implied by Eq.([T7|l is that the absolute 
value of the critical current is independent on the direction of 
the current, i.e. - with 



/+ = max 7(0) , /; = |min 7(0)1 . 



(18) 



Over the next Sections, we discuss how the barrier presented 



in Sec II B provides a controllable way to violate the relation 
- . In this section we discuss the fate of the properties 
Eq.([T7]i when Eq.([T5|l no longer holds due to the presence of 
external magnetic fields coexisting with spin-orbit coupling. 

Starting from the zero magnetic field condition we group 
the Andreev states in 7"-related states, \i) and | j), and address 
the response of the energy levels to a weak external magnetic 
field. As I j) is the time-reversed state of they have opposite 
local spin projections. 



{Sj(r)) = {j\s(r)\j) = -0>(r)|J> = -<sj(r)) , 



(19) 



with s - jcr and cr - (a-x,a-y,cr,). Let us assume, as an ex- 
ample, that the spin of the Andreev state \j} is up along the 
y-direction — then we know that |j) is a down state. Thus, if a 
Zeeman field is applied along the y-direction — provided it is 



not too strong to alTect the properties of the original states — 
the energy levels corresponding to \j} (at <p-(px) and to |j) (at 
<f) = 2n—(px) change in the same amount, \ez\, but with opposite 
signs. In general, if the total spin of the Andreev state is not 
fully aligned with the external magnetic field then the splitting 
will be a fraction of |ezl (but still proportional to the magnetic 
field). We must also keep in mind that the states change their 
spin structure as (and e) varies. 

In order to qualitatively understand how phase-shifts may 
arise due to the Zeeman energy shift, sz, we now particularize 
to the case of a transparent N-S interfaces in a long junction — 
a situation in which their Andreev levels ca n be a pproximated 
as linear functions of the phase for |e| < Ao.^^For |£| < Aq 
we write the levels as 



e/0, ez) = Aj{4> -n)-ez-{S ;,,.) , 

2 

ej(0, ez) = -Ajicp - n) - sz-{SjJ , 

2 

= -'lj{4>-7r) + sz-{S ;;,.), 



(20a) 
(20b) 
(20c) 



where Aj > is the slope of the level and (S ^,y) is the total spin 
along the y-direction of the Andreev state \q} — we later re- 
late these quantities with the parameters of our system. From 
Eq.(|20]i we see that the shift in energy, due to the presence of 
the magnetic field By, is equivalent to a shift of the supercon- 
ducting phase of the Andreev levels at zero field. 



e/0,£z) = fi/0-0),O), 
ej(0,£z) = fij(0-0f,O), 



■(Sj,y). 



(21a) 
(21b) 

(21c) 



The energy levels in Eq.pO]) and Eq.(|2T]i are single-valued 
functions in 0. However, by translating all solutions to the 
[0, 27r] region, we get many Andreev levels for each value of 
in that range. The effect of the magnetic field is shown in 
Fig|5ja) for an unrestricted phase axis while Fig|5jb) shows 
the same results folded into the [0,2;7r] region. This illustrates 
how the energy levels of a T-pak of Andreev states would 
rearrange due to the magnetic field in any long junction in the 
transparent interface limit. 

From the above discussion we see that the 7~-pairs provide 
a useful way to group the Andreev levels when discussing how 
a weak magnetic field can modify them. It is appropriate to 
make an important remark here. It is known that for a sys- 
tem with zero spin-orbit coupling but nonzero magnetic field 
the Josephson current fulfills the constrictions of Eq.([T7]l ir- 
respective of the transparency of the interfaces, the value of 
the magnetic field and the length of the junction. This fol- 
lows from a symmetry of the Hamiltonian analog to the one 
presented in Eq.([T5]l for time-reversal systems, namely,'22l 



(22) 



where U - TR„ and R„ is a 7r-rotation of the spin along 
an axis perpendicular to the one of the magnetic field. The 
TR operator, T, reverses the magnetic field direction and R,^ 
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FIG. 5. (Color online) Phase shifts due to external magnetic fields. 
The N-S interfaces are transparent, (a) Andreev levels in the long 
junction limit. At zero field, due to the symmetry under T of the 
BdG Hamiltonian, for every state at there is an associated 
state |7>, at In - (py (= -^v) that: (i) has the same energy, (ii) has op- 
posite local spins properties (see text), and (iii) has opposite velocity. 
A weak in-plane magnetic field By — if applied along a direction in 
which the total spin of the Andreev state is nonzero — generates an 
increment (decrement) of the energy of the state \ j) as shown 
with vertical arrows. The new pattern of Andreev levels (thick lines) 
are phase-shifted (horizontal arrow) from the S = levels (thin lines), 
(b) Shows the Andreev levels of (a) restricting (j) to the (0, 2n) region. 

maps the field back to its original direction. The relation holds 
even with a position dependent magnetic field, providing the 
magnetic field is coplanar — the rotation R„ is taken along the 
normal to such plane. The properties of Eq.([T7]i for the CPR 
follow straightforwardly by noticing that the two states in a 
IZ-pair (analogous to states in the 7"-pair of Eq.([T6]l) have 
opposite velocities. 

The important point is that, for Josephson junctions in 
which SOC and Zeeman fields coexists, a symmetry opera- 
tor U' as the one in Eq.p2|) does not exists in general and 
thus the CPR properties shown in Eq.([T7ji can be violated: the 
CPR can have 7(0) 0, 1{n) + and /+ + I'W It is instructive 
to discuss, for example, the case of ID spin-orbit system with 
a Zeeman field, 

•W" = -^SPx - nk,(rxf -E, + Ez,icri + Ez,i(r^ , (23) 
Zm 

where we have written the {1,2, 3)-directions in spin space for 
gaining generality (cr, are the Pauli matrices) with the SOC 
along the 1-direction. Ez,\ and £23 are the parallel and per- 
pendicular Zeeman energies, respectively. This Hamiltonian 
can be mapped by a well-known unitary transformation to a 
system without SOC but with an inhomogeneous magnetic 
field. The transformation is just the spin-dependent shift in 
momentum,^" Ua - exp (i^„;iccri). In the transformed Hamil- 
tonian, the parallel magnetic field is unaffected, whereas the 
perpendicular component rotates along the 1-axis as a func- 
tion of the position x, 

Ul^'^U,^ ^^-E„+Ez.i(Ti 

+Ez,3 [cos (2kaX) 0-3 - sin (2kax) 0-2] . (24) 

Notice also that the s-wave superconducting pairing potential 
transforms trivially under Ua-^^ This simple argument allows 
us to see that for magnetic fields that are neither parallel nor 
perpendicular to the SOC axis — i.e., Ez^i + and £23 + — a 
symmetry as in Eq.(|22|) cannot be found because the system 



is intrinsically equivalent to a system without SOC but subject 
to a non coplanar Zeeman field texture. 

For the 2D case with Rashba SOC, it has been shown in 
Ref. |26]that a symmetry tl' that assures zero Josephson cur- 
rent for = can be found only if the Zeeman interaction is 
perpendicular to the plane of the SOC.^^ For inplane fields, 
as Eq.([T7|) no longer holds, current for = and = tt is 
not forbidden by symmetry and also one may find that l^+l'^. 
However, the geometry of the sample plays a crucial role since 
not any in-plane direction is efficient to modify the Josephson 
current. For instance, due to the sample configuration (width 
going to infinity) phase-shifts in Ref. 12 are not observed for 
both directions of in-plane magnetic fields. Similarly, in the 
finite width stripe geometry we investigate,^^ (with or without 
a symmetric barrier inside) if the magnetic field is parallel to 
the stripe the current for = becomes zero. 

From the above discussion we have deepened our physi- 
cal understanding of the presence of a finite Josephson cur- 
rent at zero phase difference. Our TR-symmetric device is 
able to break the spin degeneracy of the Andreev states in 
absence of magnetic fields. Due to the sample-specific ge- 
ometry, a given in-plane direction is privileged in the spin of 
the Andreev states. More specifically, with a single transverse 
mode open inside the barrier the two Andreev states with pos- 
itive velocity (|A) and \B) states in Figs |3] and |4]) can have the 
same sign of total magnetization along the direction in which 
the SOC device polarizes. From the concepts that led us to 
Eq.([2T| it follows that a nonzero current for zero phase dif- 
ference is expected if a weak magnetic field By is applied 
as 0^ + -0|. On the other hand, current for zero phase dif- 
ference is not expected for inplane magnetic fields along the 
x-direction because for those states the integrated magnetiza- 
tion along x-direction is zero, i.e., {S a,xI - {S b,x} = 0. This is 
in agreement with our numerical simulations. 



3. Andreev levels and phase-shifts in the weak magnetic field limit 

In order to analyze our numerical results quantitatively, we 
now relate the two essential features of the Andreev states in 
transparent junctions, namely their phase dependence (slope) 
and their energies as a function of a magnetic field (in lin- 
ear response), with the parameters of our system. The ener- 
gies of Andreev states for a homogeneous and long (Lat » 
^iv = hvp/nA) one-dimensional S-N-S junction are given by 
e"''' = fivp/{2Lff)[27t(n + l/2) + (l)] where vp is the Fermi veloc- 
ity of the normal material, is the length of the junction and 
the plus and minus signs correspond to the two signs of the 
excitation velocity. This expression can be generalized to in- 
clude both the effect of external field and corrections of order 
^at/L/v. Following the procedure of Ref.l25]we obtain 

4'^ = cru^ffBy + [2n{n + 1 /2) ± 0] , (25) 

1{Lm + TTtM) 
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here cr is the spin of the Andreev state and yUefj is the effective 
magnetic moment 



1 + 



gs_ 
Sn 



1 



1 + ^ , 



(26) 



We assume that the magnetic field is in the y-direction — since 
it is contained in the plane of the 2DEG no diamagnetic effects 
are induced. The energies of the Andreev states are then sim- 
ply shifted by the field. In Eq.(p5|) n is an integer and is 
restricted to the [0, 2;r] region. These energy levels are equiv- 
alent to the ones presented in Eq.(j20| after folding them into 
the [0,2;7r] interval, as done in Fig.lsfb). 

This analytical result for a ID system without SOC can be 
used to qualitatively analyse the effects induced by the pres- 
ence of a barrier in a normal region when a + ^. We re- 
sort to the WKB approximation and proceed in the following 
way: for an adiabatic QPC in the absence of external mag- 
netic fields and for an electron (and the reflected hole) in the 
channel / we define a 'local' wavevector, 



(27) 



These functions have to be obtained from the dispersion rela- 
tions at each position x^, by assuming that the potential for all 
X is equal to Y{xa,y)- As shown in Sec |llB| due to the SOC 
the different transverse modes with dispersion relations given 
by Eq.(|7]) become mixed and the form of Ff''(x,s,fi) be- 
come nontrivial — in the absence of SOC, different transverse 
modes do not mix and the functions are just Ff''(x,e,iu) - 



^^(yu ± e) - Vi(x), with y,(jc) the effective ID potential that 
accounts for the shift of the bottom of the /-channel band. 

The orbital phase accumulated by the electron and the re- 
flected hole in the junction is given by 



(k';(x) - k'l(x)) dx 



2Lf^6iE 



(28) 



where we have parametrized the integral by a single dimen- 
sionless number 5,. With this parametrization we find for the 
energy of the Andreev states 



[2n(n + l/2): 



(29) 



Within this picture, the effect of the barrier only introduces an 
effective length L^J, > L^i for each channel — the value of vp 
is taken far away from the barrier potential. There is also a 
small correction to the effective magnetic moment of channel 
/ that is given by Eq. ( |26] l with replaced by the effective 
length Lf^Si. 

We found, by inspecting the Andreev levels's slopes of 
the exact numeric results, that the parameters 5, are only 
weakly dependent on the SOC (see below). As the spin 
is no longer a good quantum number, in the above expres- 
sion the spin index should be replaced by a new quantum 
number that characterizes the two branches of each of the 
i = 1,2, ...,A^of transmitted channels. In order to include 



the effect of an external magnetic field in our simple WKB 
picture, we profit from the properties of the Andreev states 
introduced in the previous section. We group the Andreev 
states in 7~-partners, \ j} and \ j}, given in Eq.([T6| [see Fig.js), 
where (;, J) = (1, 1), (2, 2), {2Noc, 2N^,), and ; (J) labels the 
states with positive (negative) slope. Therefore, in the linear 
response regime, a magnetic field By leads to the phase shifted 
energy levels of Eq.(|20]l. In fact, the WKB version of Eq.pO]) 
can be written as 



e/0) = rijU^fiBy + 



2(Ln6 j + ti^n) 



2{Lti6j + n^N) 



[2n{n + \l2) + (l)], (30a) 
[2;r(«H-3/2)-0], (30b) 



where the parameter \rij\ < 1 quantifies how different from the 
spin degenerate case the system behaves. The link between 
the parameters and the phase-shift 0^ defined in Eq.(|20]l can 
be written as 



A, 



\T]jfXeSBy\ = 



ez^{S j,r) 



,Z r, 2(i'A'5j + 71^ n) 
Pj - '?./^eff^v 



(31a) 
(31b) 
(31c) 



The factor n^^ in the denominator can be neglected in the 
long junction limit — this is the case of the results presented 
in Fig.|6]that we discuss below. Introducing the linear disper- 
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FIG. 6. (Color online) losephson current at a fixed finite Zeeman 
energy, Ez = /JeffSv, and a fixed barrier profile — Vg is set to allow 
transmission of a single transverse mode — as Ljy is enlarged. For 
convenience we plot it as a function of the eff'ective length, L^g = 6 jLn 
(see Eq.l |31c[ l for fjv ^ ^jv)- Panels (a) and (b) correspond to o- = 0, 
whereas panels (c) and (d) to a ^ 0. Qualitatively similar results 
are obtained if L/v is fixed and Ez is varied. For a = the normal 
device behaves as a ferromagnetic polarizer; the current at = is 
zero, while changing the length produce 0-;r junction transitions. On 
the other hand, for the case shown for a t 0, the QPC behaves as 
an efficient polarizer in the normal state. The magnetic field By is 
applied along the SOC polarization axis. Note that the current at 
= is nonzero, and that the phase-shift changes with the length of 
the junction. 
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sions of Andreev levels in Eqs.([T2|i and ( [T3] l, the zero mag- 
netic field and zero temperature contribution to the Josephson 
current of each ID mode (a set of j an d /' pa rtners) is propor- 
tional to the saw-tooth function, Fsti^)^^^ 



2eAj 



2eAj (f, 



if < < TT 



^(^-2) ifn<4,<2n ' 



(32) 



which is a 27r-periodic zero-mean function and we have de- 
fined it for (f> e [0, 2ti\. For finite the total current includes 
the contributions of the Andreev bound states (|e| < Aq) and 
of the continuum states (|e| > Ao)p2lFor a perfect ID channel 
the CPR has an abrupt change ?A (p-n, however, in the sim- 
ulation such jump is slightly smoothed due to the presence of 
the barrier (see Figj3|c)). 

We now discuss the effect of the magnetic field on the 
Josephson current. The gate voltage Vg is set to allow the 
transmission of a single channel, i.e., Noc - 1- Thus, there are 
two distinct Andreev levels, {j, j) - (1, 1), (2, 2), with positive 
and negative slope, respectively. In the simulation there is a 
weak in-plane magnetic field (i.e., Ez <K Aq) and we vary the 
length of the junction while the barrier dimensions are kept 
fixed. The quantities SjLn change and so the phase shifts in 
Eq.pia[). The current becomes 



if) + -^F,M - , 



if)- 



(33) 



We show the CPRs for a system without SOC in Fig. |6|a) and 
|6|b). Note that there is no current for (p-Q because the states 
i - \ and i - 2 are the spin up and down states along the y- 
direction, respectively, and so 0f = -(/>f (because 772 - -rji - 1). 



For 



^f = -0f- 
functions is at 



(2n + \)n the discontinuity of both saw-tooth 
p - 0. This can be understood as a transition 
from a to a ;7r -junction. 

Figure |6jc) and |6|d) show the CPR for a system with po- 
larizing properties due to the presence of strong SOC. The 
current for (/> = is nonzero as both saw-tooth components 
move in the same direction. This can be understood as fol- 
lows: as the mode (1, -h) dominates the Andreev state |A) [see 
Fig. |3je) and Fig. QA)] we have 771 = t/a ~ 1- On the other 
hand, since the region of the barrier is small compared to the 
total length of the junction, the (2, +) mode dominates away 
from the barrier in the Andreev state \B) [see Fig. [3|f) and 
Fig. |4|B)], and we also have 772 = ?7b ~ 1- It then follows that 
02 justifying the numerical results of Fig. pic). Using 
Eqs.((33]l, ( |32l ) and ([30| we see that the longer the junction 
the bigger the phase-shift (oc L^) while the maximum current 
gets reduced (because Aj oc L^^^). Similarly, if is fixed 
both the phase-shift and 7(0 = 0) grow linearly with B, while 
the critical cuiTent does not change. 



4. Phase shift of Andreev levels, dependence with a 



are fixed. As expected, for a-0 [Fig. |7]^Lc)] the arrangement 
of Andreev levels is symmetric with respect io (p-n. We call 
(j' y) = (1 > 1) \ih i) - (2, 2)] to the two energy levels that cross 
£ = at > TT [at < ;7r]. Notice that as a is increased the 
position of the crossing for (2, 2) at e = shifts in phase to the 
right whereas for levels (1,1) the crossing position remains 
unchanged. This is interpreted using the WKB picture . We 
associate levels (1,1) to states in which electrons and holes 
do not go through avoided crossings in the transport through 
the QPC. Therefore, in these states the spin remains along the 
y-direction: i.e., 771 =s -i-l, for the Andreev state associated to 
I A) [see Fig.|4jA)] and its T -partner. 

On the other hand, the Andreev states corresponding to lev- 
els (2,2) have their spin not fully aligned to the y-direction 
as electrons and holes are affected by the avoided crossings 
for nonzero a. In the strong SOC limit we associate them to 
the Andreev state |B) [see Fig. QB)] and its 7"-partner. As 
the probability to follow the avoided crossing increases with 
a [see for example the 2-channels model in Eq.(|9|] the wave- 
function in the region away from the barrier is changed from 
being dominated by the first transverse mode (at a a; holding 
772 = -1) to being dominated by the second transverse mode 
with opposite spin along the y-direction (at very strong SOC 
such that 772 ~ 1). 

For both types of states we take the slopes of the Andreev 
levels Aj and from Eq. pia[ ) compute the values 6j. Notice 
that we have sorted out the 7~-partners (here j - 1 , 2) by in- 
creasing i as Aj decreases, i.e., from the fastest to the slowest 
Andreev state. Then we introduce the observed shifts 0f (see 
Fig.|7|Lc)) in Eq.lpO]) to compute the parameters 77^. The re- 
sults are presented in Table[I]and Fig.[7|Lb) where we plot the 
values T]j. 

In Fig. |7|Lb) we see that, as expected, 771 is constant as 
a function of a. We show that 772 is well approximated by 
2Pav - 1, with Pav the average polarization of the QPC in 
the E € (-A(), Aq) energy window. This quantity is a crude 
estimate for the spin projection of the Andreev state that 
suffers the avoided crossing for the case in which the total 
length of the junction is much larger than the zone of the bar- 
rier in which the spin gets mixed. In this limit one can as- 
sume that the spin of the Andreev state is dominated by their 
spin properties away from the barrier, these being the ones 
which in the normal device dictate the value of the polariza- 
tion. More specifically, we assume that 772 ~ 5 2, with Si ex- 
tracted from the polarization of the normal device. We assume 
^'av = (1 + S2)/2: the factor 2 in the denominator accounts for 
the total number of transmitted modes within the first conduc- 
tance plateau and the numerator is the sum of 1 (i.e., full spin 
up due to the electron never leaving the (1, +) mode) with the 
contribution S 2 given by the specific superposition in modes 
(1,-) and (2,-1-) arriving to lead R due to the SOC-induced 



avoided crossing (see 2-channels model in Sec 11 B 1 



In Figure |7|L) we show the Andreev levels and the CPRs 
for several values of a corresponding to the weak, intermedi- 
ate and the strong SOC regime. The number of open channels, 
Noc, is one and both the junction length and the magnetic field 



5. Asymmetric CPRs 

So far we have shown results for one open channel in the 
QPC. In the weak-B, limit we have shown that a current ap- 
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FIG. 7. (Color online) Anomalous current, /(0 = 0) # 0, phase-asymmetric CPRs and current-asymmetric CPRs, due to in-plane magnetic 
field (parallel to the SOC polarization axis). The barrier is the quantum point contact, QPCi, described by Eq.( |Al[ l taking (z; Wi,;Lh) = 
(30; 100; 250)nm. The junction length is Ljv = 1 .2;um and the S-N interfaces are transparent. The Left (L) panels show CPRs (La) and Andreev 
levels (Lc) for a fixed magnetic field, Ez = gi-iBBy = = 150/jeV, and difi"erent values of the Rashba strength, a. The gate voltage, Vg =46.5meV, 
sets the number of open transverse channels, A'^^, in 1. The larger a the bigger is the current I((p = 0) and the less distorted is the Josephson 
current from the one at fi^ = 0: it looks as a phase-shifted version of the CPR at zero field (see (Ra) and (Rc)). (Lc) Deduced relative magnetic 
moments, and fi2, of the two Andreev bound states states (see text) as a function of a. We show the simplest estimation for ji2, (2Par - 1), 
with P„y the average polarization induced by the device in the energy range of the superconducting gap. The Right (R) panels show how the 
opening of more than one transverse mode in the QPC can lead to asymmetric CPRs when By is nonzero. We fix Q- = 40meVnm and show the 
Andreev levels, for (Vg,N„,) = (46.5meV, 1) [ (39.4meV, 2)] in panel (Ra) [(Rd)] for fi,, = and in panel (Rb) [(Re)] for B,. tO: Ez= ISOyueV. 
In panel (Rf) the CPR for By 1 and A'^^. = 2 does show current-asymmetry, i.e., it fulfills Eq. ^Ijj. On the other hand for the case of Noc = 1 
shown in panel (Rc) the CPR gets phase-shifted and distorted for BytO but no strong current-asymmetry is observed. 



TABLE I. Dependence with a of parameters 6j and r]j, j=l,2. The 
JJ has the quantum point contact QPC] , the CPRs and Andreev levels 
are those shown in Fig.|7](left panels). 



TABLE II. Values of 6j and j for two open channels in the QPCs for 
different junction lengths when Q' = 40meVnm. 
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-0.90 


-0.77 


-0.44 


-0.03 


0.34 


0.59 


0.74 



pears for zero phase difference. This is due to different shifts 
suffered by each of the two 7~-partners of Andreev states. As 



we discuss in Sec III A 2 here the current-phase relation does 
not need to be symmetric and thus the critical current can de- 
pend on the direction of the current flow, i.e., I^. We 
find that the latter effect becomes significant when the barrier 
is tuned to allow the transmission of a few (more than one) 
open channels. This is shown in Fig. |7|Rf) for Noc = 2. For 
convenience we number the 7"-couples y = 1,2,..., 2Noc from 
faster to slower Andreev states (decreasing slope of their de- 
pendance on <p). In the case of Fig. |7|Rd) and Fig. |7jRe) we 
show how four of those couples (Noc - 2) rearrange due to a fi- 
nite magnetic-field applied along the direction of polarization. 
In Table [ll| we present the parameters 6j and rjj obtained from 
the simulations. 

As expected the fastest 7~-pair (7=1) has 771 * 1 due to the 
fact that they are associated with Andreev states in the out- 
ermost branches of (1,-1-) and (1, -), i.e., those unaffected by 
avoided crossings. The remaining three 7~-pairs are affected 
by the avoided crossings in the transport through the barrier 



Vg [meV] LA.[nm] 


Si rj, 


62 ri2 


S3 m 


^4 '74 


QPCi 600 
39.4 900 
1200 


1.39 0.95 
1.27 0.97 
1.21 0.99 


1.49 0.50 
1.35 0.55 
1.27 0.58 


1.63 -0.31 
1.44 -0.41 
1.33 -0.45 


1.82 0.18 
1.58 0.27 
1.47 0.36 


QPC2 600 
46.3 900 
1200 


1.53 0.97 
1.40 0.97 
1.31 0.99 


1.60 -0.84 
1.45 -0.85 
1.36 -0.86 


1.70 0.87 
1.52 0.89 
1.41 0.90 


1.90 0.72 
1.66 0.77 
1.52 0.79 



Each of them also sees a different effective length (djLf^) in 
the passage through the barrier This combination of differ- 
ent slopes and different spin properties (rij) translates, given 
a finite magnetic field along the y-direction, into an asym- 
metric CPR. The asymmetry is difficult to detect in the same 
device for the Noc = 1 case (at a higher gate voltage). In 
Fig. |7|Ra) and Fig. |7]^Rb) we show that the rearrangement of 
Andreev levels for Noc = 1 is symmetrical (with respect to a 
new phase value different from n) except for the small differ- 
ence in slopes between the two 7~-pairs; the latter difference 
does not produce a significant asymmetry in the total CPR 
shown in Fig.jTj^Rc). 

The significant difference of the two critical currents, |/^ - 
Icl/i^c )' ^Iso gets reduced in the opposite limit, when 
Noc s> 1. In this case the full current is the sum of the contri- 
butions of the 2Noc 7~-pairs and balanced spin-properties are 
dominant since a smaller fraction of the transmitted channels 
undergo an avoided crossing transition. This can be under- 
stood, for example, by taking the height of the barrier half 
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of what is considered in Fig. [T] Then, the electrons in chan- 
nels (1 , -) traveling from left to right would not be affected by 
the avoided crossings with mode (2, +) and thus both modes 
(1,0") would be transmitted. Assuming that the third trans- 
verse mode can now transmit, only the SOC mixing between 
modes (2, cr) and modes (3, o") contributes to spin filtering ef- 
fect. This is the reason why the polarization, and thus the 
anomalous Josephson effects we study here, is smaller when 
many channels are allowed to transmit through the barrier. 



in Fig. Re)) 



Asymmetry = 



(36) 



The underlying resonant tunneling physics produce an abrupt 
dependence with the gate voltage of the quantities of inter- 
est. This is different from the idealized S-N junctions case 



of Sec III A for which the results do not change with Vg pro- 



Anomalous Josephson efTects enhanced by 
barriers with SOC: non-ideal interfaces 



In the previous section the analysis was deliberately made 
to elucidate the effects introduced by the barrier. Therefore 
we have idealized the S-N junctions in order to minimize nor- 
mal reflections at the S-N interface. In those simulations the 
matching between 2DEG and superconductors was very good 
despite the fact that the latter have or = and a smaller g-factor 
Her e we ch ange the superconductor system, as described in 
Sec 



II A 



22| , 



with the goal of incorporating the mismatch of 
Fermi velocities between the 2DEG and the S contacts to sim- 
ulate more realistic experimental devices.^' In the following 
we show that anomalous effects in the CPR due to the polar- 
izing effects are present both for realistic junctions based on 
quantum point contacts and for wide barriers. 



1. Quantum point contacts 



The drastic effect of more realistic 2DEG-contact inter- 



faces was discussed in Sec II B 3 for the transport properties 
of QPCs when the leads were taken to be metallic, i.e. Aq - 0. 
As shown in Fig. [2] — for a ratio between Fermi velocities, 
Vp^'^/Vp'' si2.25 — the conductance through the QPC shows a 
resonant-like behavior as a function of gate-voltage. In partic- 
ular we have shown that, due to quantum interference effects, 
the polarization peaks can be larger than the polarization de- 
veloped for the case where the 2DEG-contacts interfaces do 
not generate backscattering. 

Figure [8] show results derived from the current-phase rela- 
tion in a Josephson junction containing a QPC and realistic 
2DEG-S interfaces. Due to the magnetic field By, the CPR 
for the system has interesting properties that can be tuned by 
changing the gate voltage which controls the barrier height. 
For convenience we define the average critical current (see 
Fig.la)) 



(34) 



the anomalous current (or zero phase difference current, see 
Fig.gb)) 

/„ = 7(0 = 0), (35) 
and the asymmetry of the CPR (see the dependence with Ez 
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FIG. 8. Results for a Josephson junction, S-QPC-S, based in QPC2 
with realistic 2DEG-S junctions: an In-based 2DEG and Nb super- 
conductors are simulated (parameters given in Fig. |2jb)). We take 
a = 5meVnm, Ln = 1 .2yum; in (a) and (b) for the JJ results the Zeeman 
energy in the 2DEG is Ez = -150yueV. We plot the dependence with 
the gate-voltage, Vg, of: (a) the average critical current, /,., and (b) the 
anomalous Josephson current !„. In (a) [(b)] we include the conduc- 
tance G [the polarization P] as a function of V^, for the same device 
between normal contacts (dashed line, 2DEG-QPC-2DEG junction) 
such that the leads do not generate backscattering, and the leads are 
metallic (solid black line, M-QPC-M junction) we take the parame- 
ters of Nb with Ao = 0. We observe correlation between 1^ {I„) for the 
S-QPC-S and G (P) for the M-QPC-M. (c) We fix Vg at a peak for 
/„ and present, as a function of Ez, L and the critical current asym- 
metry For such value of the apartment from the I* = 1^ 
situation is strong and controllable via magnetic field (and gate volt- 
age as shown in Ref. 22). (d) We plot the full current-phase relations 
for some situations labeled in (c). 
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viding the number of open channels in the QPC remains the 
same. The dependence with Vg shows resonant peaks in both 
Ic and /„. Fig. Sjd) presents some examples of CPRs that show 
a striking asymmetry. In these cases, the fact that the CPRs 
fulfill the zero area condition of Eq.([T4]) is not obvious. 

Indeed, because of the strong normal reflection at the in- 
terfaces, the system can be interpreted as a S-quantum dot-S 
Josephson junction and thus the resonant like-shapes in the 
critical current are expected.!^ This is the so-called resonant 
tunneling supercurrent transistor regime for samples in which 
the charging energy of th e "do t" is small and thus Coulomb 
blockade is not observed.l^SH jn (jjg absence of a magnetic 
field the states arising in the quantum dot appear as time re- 
versal pairs. As for the case of Sec III A for JJs based on 



normal quantum dots with spin-orbit coupling, no asymmet- 
ric CPRs or current for zero phase diff'erence are expected at 
zero field.l^ 

A correlation between the normal transport and the super- 
conducting transport is observed in the shapes of Ic and G 
in Fig. [8|a). This is understood generically, for Josephson 
transistor devices, because the critical current is dominated by 
the contribution of the highest quasiparticle level (i.e., clos- 
est to e = 0).^^ Then the critical current is maximized at the 
same values of Vg that maximize G in the M-QPC-M device: 
constructive interference favors transmission of quasiparticles 
and Andreev bound states acquires strong dependence with 
the phase for e a; 0. Moreover, in Fig. [HJb) we find a strong 
correlation between P in the normal device and the anomalous 
current in the associated Josephson junction. 

Here the fact that the transport is conditioned by normal 
backscattering at the interfaces and interference effects does 



not invalidate the symmetry arguments discussed in Sec 111 A 



for the existence of Andreev states 7"-partners at zero mag- 
netic field. Indeed, the Andreev states at a resonant condi- 
tion are 7~-partners and it is natural to expect that their spins 
(see Eq.([T9]l) are correlated with the polarization found in the 
corresponding normal junction. At the resonant value of Vg 
their associated Andreev levels cross each other with oppo- 
site slopes at (e, <p) - {0,7t).^'^ In a zero SOC system Andreev 
states are spin degenerate and thus there are two 7"-pairs res- 
onating at the same Vg. Here on the other hand, as pointed out 



in Sec II B the two scattering states contributing to a given 



transverse mode have different velocities. In the Fabry-Perot 
geometry this leads to resonances at different values of Vg. 
In summary, the multichannel physics, the symmetric sample 
geometry, and the SOC make a particular spin projection, up 
along the y-direction in this case, to be favored in the normal 
transport from left to right (Pl^r > 0). The opposite pro- 
jection is favored in the right to left normal transport. The 
Fabry-Perot effects enhances the polarizing effect for specific 
values of Vg. In the associated Josephson junction, when a 
weak By breaks the spin degeneracy, this leads to a nontrivial 
rearrangements of the Andreev levels resulting in the nonzero 
la and/or asymmetric CPRs, i^I^. 

It is worth mentioning that the biggest asymmetry in the 
CPR occurs due to the magnetic field for the cases that at zero 
field present a few off-resonance 7"-pairs (i.e., pairs whose 
energies are below the Fermi energy for any value of the phase 



difference) and a single resonant 7~-pair (its energy reaches 
e = for = tt). In this situation, the dispersion of the res- 
onant Andreev state is close to the one corresponding to the 
ballistic case while the off-resonant ones are more sinusoidal- 
like, typical of a low transmission channel (results not shown). 
Therefore, the shift induced by the B field are very different in 
each case, leading to the observed asymmetry when the cor- 
responding contributions to the CPR are added up. This be- 
comes strongly dependent on Vg allowing, at a fixed B, to use 
slight changes in the gate voltage as a control parameter^^ The 
condition is generically observed as long as few open 
channels are open in the QPC and the modes affected by SOC 
mixing contribute a significative fraction of the total current. 
We observe that the asymmetry grows with the magnetic field 
as shown in Fig. |8jc) up to a value in which it jumps from 
> 7^7 to < 7^7. This dependence is also observed at fixed 
By as a function of the junction length — in such a case one 
has to change Vg with in order to follow the resonant con- 
dition. The qualitative interchangeability of and By, as in 
S-F-S junctions, follows because both quantities produce an 
enhancement of the underlying phase-shifts of Andreev levels 
(see Sec |IIIA3| l. 



2. Wide barriers 



In Sec II B we have presented thepolarizing properties of 
QPCs.^^ The WKB picture, in Fig. fTl shows therefore trans- 
verse modes well separated in energy by virtue of the transver- 
sal quantization induced by the lateral constriction of the 
QPC. The polarizing properties, however, are also expected 
for very wide barriers. The physical mechanism is concep- 
tually the same, i.e., a smaller transverse modes quantization 
energy does not avoid the SOC-induced mode mixing leading 
to polarization in the transport through an adiabatic barrier 
This nontrivial spin transport shares the same time-reversal 
properties as in the discussed QPC. Nonzero polarization has 
been derived by Silvestrov and Mishchenko in Ref. 29 even 
for the case of a very wide system with transversal periodic 
conditions. As the barrier (B) becomes smaller the polariza- 
tion decreases in agreement with the observations in QPCs. 
Here we present results for the polarization in a wide barrier 
and show that these devices also lead to the same CPR phe- 
nomena in Josephson junctions. 

Figure [9] shows results for the wide barrier potential 
Vv,ide{x,y) given in Appendix [A] The barrier has no depen- 
dence with y and thus lateral confinement is determined by 
the stripe width Wy chosen in the simulation. We take Wy = 
630nm, which means 51 open channels in the stripe (V^ - 0) 
for our choice of the Fermi energy. We denote with Z = to 
the case of a 2DEG-B-2DEG junction, i.e., the leads are es- 
sentially the same material as the central region except from 
the reduced g-factor and the zero spin-orbit coupling-recall 
that Z is the amplitude of the barriers at the S-N junctions 
in the BTK approach.-22i Then we study the case in which the 
leads are metallic, i.e., a M-B-M junction. In the latter case we 
use the parameter A presented in Eq.(|6]), to reduce the trans- 
parency of the junction away from the ideal condition (A = 1). 
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In Fig.|9|a) we shown the normal transport properties of the 
barrier, the conductance G and the polarization P, that we use 
to characterize it. In the Z - Q case, despite the fact that the 
conductance plateaux become undefined, the polarization is 
nonzero. In the case of metallic leads with high transparency, 
A-1, the G and P acquire some structure but without show- 
ing marked resonances as in the QPC case. This can be under- 
stood from the fact that many channels with similar velocity 
are available and thus neighboring resonances overlap. The 
total contribution is similar to the Z-Q case except for a re- 
duction of the conductance. When we take A = 0.3, the trans- 
parency of the junction is decreased and both the conductance 
and polarization are affected. In particular the polarization 
shows peaks as a function of Vg. 

We transform all these cases into Josephson junctions by 
taking Aq = 1 .5me V. The current-phase relations in the pres- 
ence of a finite magnetic field along the y-direction is shown 
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FIG. 9. Results for a junction with potential V„iic{x,y) with = 
360nm, i.e., a plain barrier (B) without a lateral constriction. We 
take a = 20meVnm and the 2DEG and superconductor parameters 
given in Fig.[2]and Fig.jS] (a) Spin polarizing properties of the barrier 
when the leads are normal, we show the dependence on gate voltage 
Vg for the conductance G and the polarization P. Case Z = is for 
a 2DEG-B-2DEG junction, i.e., the interfaces are ideal as the leads 
are the same In-based 2DEG. The other two cases are for the M- 
B-M junction, the leads are metallic Nb and a there is mismatch of 
Fermi velocities and effective masses; moreover, away from A = 1, 
see Eq.l|6j, the transparency of the junction is decreased, (b) The 
normal systems of (a) are turned into Josephson junctions due to a 
superconducting pairing Ao = 1.5meV in the contacts and we also 
include a Zeeman field along the y-axis with energy Ez = -150/ieV. 
Current-phase relations showing both current for = and/or CPR 
asymmetry tl^, as those appearing in QPCs, can also be obtained 
for wide barriers. 



in Fig.|9|b). We can easily find, by tuning Vg, situations with 
la + and asymmetric CPRs. As it is the case for QPCs, the 
asymmetry of the CPRs becomes larger for the realistic junc- 
tions having normal scattering at S-N interfaces. The phenom- 
ena is smaller the more modes are allowed to transmit in the 
device as it is also the case for the polarization of the associ- 
ated normal device. 



IV. CONCLUSIONS 

We have described how the spin-polarizing properties of 
barriers with spin-orbit interaction affect the spin of the An- 
dreev states. First we revisited the normal device polaiiza- 
tion mechanism highlighting the multichannel nature of the 
transport while stressing that the polarization of the transmit- 
ted electrons is due to the mode mixing induced by the SOC 
term proportional to the transversal linear momentum. As the 
device is time-reversal-symmetric strong restrictions are im- 
posed on the spin-resolved transmissions in opposite direc- 
tions. At zero magnetic field the properties deiived from these 
restrictions do not conspire against the formation of Andreev 
bound states, this is opposite to the case of ferromagnetic po- 
larizers. 

At zero magnetic field a symmetry of the BdG equation,'^^ 
including the SOC, links Andreev states corresponding to 
phase differences between superconductors -1-0^^ and -0^, re- 
spectively, that have opposite spin and velocity. The sym- 
metry argument is valid for the junction irrespective of the 
transparency of the S-N interfaces. Each of these states, that 
would appear spin-degenerate in the absence of SOC, are non- 
degenerate in its presence. More remarkable, the polarizing 
barrier makes the two Andreev states traveling in the same di- 
rection for (with slightly different energy due to SOC) to 
share similar spin properties. Because of the mentioned sym- 
metry there are two Andreev states for -<p,^ that travel in the 
opposite direction having opposite spin than those for +<px, 
so the system presents a spin chirality. In our system we have 
shown that the introduction of a weak magnetic field along the 
direction in which the Andreev states are spin chiral is enough 
to trigger anomalous Josephson effects such as current at zero 
phase difference and . The chirality condition here is 

built into the device properties by the presence of the barrier 
and thus it can be switched on and off with gate voltage. The 
condition is similar to the one required for nonzero anoma- 
lous cuiTent in the spin-orbit coupled quantum dot Josephson 
junction of Ref. 18 As this spin chirality is intrinsic to other 
helical systems, we expect that a similar effect could occur 
in those cases too. That is, a Zeeman induced I„ would also 
appear if the B field is aligned with the helical dir ection in 
Josephson junction based on s-wave superconductors-^^l^^ 

In order to see how the spin-polarizing properties of the 
baiTier determine the Andreev states and affect the Josephson 
effect, we have idealized the S-N inteifaces to minimize nor- 
mal scattering. The numerical results agree with a simple pic- 
ture for the transport through the barrier, based on the analysis 
of the contribution of two almost degenerated Andreev states 
with the same velocity and spin for a given phase difference. 
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In this limit a direct translation of Zeeman energy shifts into 
phase-shifts of Andreev levels can be made because of the 
linearity of the phase dispersion of the Andreev levels in the 
long junction regime. This, combined with the spin chirality 
of Andreev states, leads to 4 9^0. We have studied the depen- 
dence of the phase-shifts with junction length, magnetic field 
and SOC strength. We also showed that asymmetric CPRs are 
possible when the barrier allows more than one open channel. 

We also presented results for QPC based junctions with re- 
alistic S-N interfaces. The asymmetry of the CPRs becomes 
more marked than in the transparent junction case. This is 
related to the Fabry-Perot like transport that produces peaks 
in the conductance and polarization of the associated normal 
device (taking Aq = in the leads). 

The phenomena described here are to be expected in bar- 
riers as long as the device behaves as a polarizer — due to the 
mode mixing induced by the SOC. Then, few open channels 
are preferred. We have shown that anomalous effects are also 
observable in wide baiTiers even in cases in which the trans- 
parency of the S-N interfaces are not high. Therefore for 
smooth baiTiers with few open channels in quantum wires the 
mentioned effects should be observed as long as the spin-orbit 
coupling is 2-dimensional, i.e., it must have terms both pro- 
portional to the longitudinal momentum and to the transversal 
momentum. As for quantum wires constriction-type barriers 
are difficult to implement, the latter is an important result: 
apart from the practical possibilities arising from the super- 
cuiTent rectifier and the phase shifted CPR, the experimen- 
tal observation of the effects would provide a signature of the 
presence of a SOC term proportional to the transverse momen- 
tum. We believe this is valuable information as SOC is one 
of the fundamental ingredients required in quantum wires for 
accessing the topological superconducting phase, and knowl- 
edge of all existing SOC terms is desirable. 



Appendix A 

1. Potential Barrier 

To describe the quantum point contact we use the realistic 
potential given in Ref. 56 for a split-gate defined QPC, 



VQPc(x,y) 



VgT(x - Xc,y - yc) 
T(0,0) 



(Al) 



2z ' 2z 



2z 2z 



+f{'^±zJ±. '^b-2y \ J2x + Lb Wb-2y\ 



2z ' 2z 



2z ' 2z r 



7T UV 

j(u,v) arctan u - arctan v + arctan , 

2 Vl + m2 -H v2 

where z is the distance between the 2DEG and the plane of 
the gates, Li, is the length of the gate, Wb is the width of the 
constriction, and Vg is the value of the potential at the saddle 
point located at (Xcjc). 

In the case of barriers without the lateral constriction, we 
use a modified version of the potential given in Ref. l57l 



V„ide('>C,3') 



[l+cos(^^^^)] fox -Lb<x-Xc<Lb 







otherwise 



(A2) 

where 2Lb is the barrier length and at the center of the barrier, 
Xc, the barrier achieves its maximum value V„. 



2. Tight binding model 

Following standard finite-difference procedures for a square 
lattice with spacing oq, we write the total Hamiltonian as 



H -Hz+Hjsi+Hi+Hr+Hi^ i+Hn,r ■ 



(A3) 
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For convenience we have used the same notation as the con- 
tinuos Hamiltonians given in Sec jll A| and we have introduced 
the Hamiltonians H^y connecting the normal region and the 
lead y - {L,R}. The Hamiltonian Q for the normal region 
becomes. 



reN,a- 



reN,cr 



(A4) 



where ^ (Cr.cr) refers to the annihilation (creation) operator 
for an electron in the 2DEG at position r with spin cr = cr = 
{t,i). The hopping is = fP'/2m*a^. The summations are 
taken over all lattice sites in the normal region, N; terms with 
operators Cr,o- and cl^j. with r outside such region are taken 
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to be zero. The Rashba spin-orbit coupling parameters are 
related to the Rashba strength in the continuos model as 



a 
2ao ' 




a 

2ao 


a 




. a 


2flo ' 




2ao 



(A5) 



Notice that time-reversal symmetry is intrinsic to the SOC 
coupling as we have Aj^a-s- - - (^>,(5-(7-) with j - {x, y). 
The lattice version of the Hamiltonian of the leads ([TJ is 



r,cr^r,(T ' 



reS„,cr 



Ao X 



d ^d , + h.c. 



, (A6) 



where y - {R, L} for the right and left lead, respectively, the 
summations are taken over sites belonging to the regions of 
those leads Ry and d^ ^ (dr a-) refers to the annihilation (cre- 
ation) operator for an electron in the lead y at position r with 
spin cr^ - (T = {1,1}- We use the parameters ts and Es for 
simulating the different types of leads. 

When the goal is to minimize the scattering at the leads we 
use Es =0 and ts - tf^ (and A = 1, see Eq.(|6])). On the other 
hand, for simulating realistic superconducting leads (or metal- 
lic ones if Aq = 0), we change both Es and ts assuring that the 
lead has many open channels — with the appropiated velocity 
(dispersion) — in a large energy window around the Fermi en- 
ergy. The size of the latter is not important in the metallic 
phase, as the low temperature transport is dominated by the 
quantum transmission at the Fermi energy. In the supercon- 
ducting case, it is enough to consider an energy window of a 
few Ao, as the Josephson current in S-N-S junctions is medi- 
ated by Andreev reflection which decays exponentially as a 
function of energy for |e| > Ao.'^^The contribution to the CPR 
of states at energies outside this window is negligible. 

In writing the Hamiltonians Hy we have chosen a gauge 
where the superconducting phases ((pi, <pr) - (0, -<p) and is 
accumulated only in the tunneling Hamiltonian between the 
central 2DEG and the right superconducting lead.^^ We group 
the lattice sites according to the coordinate x. The set Xr {X[}, 
having the sites at x = xr (x - xl), is the last (first) vertical layer 
in the normal region. The intermediate Hamiltonians are 



N,y 



V f^e'-^'c"'" d 



V ^e^'^'c^ d - +hc 



(A7) 



with parameters (s^, A!^^, A^^,-&l) - (-1, A*^ A*^^^,0), and 
isR,A'^^,A[^,-d-r) = (l,/l;c,ti,/lx,iT'-^/2)- Notice that the vai-i- 
ation of the Rashba strength — as a is taken to be zero in the 
leads and non-zero in the central region — introduces a SOC 
contribution to the tunneling Hamiltonian at the interfaces. 
This follows as a result of the correct discretization of the SOC 



term ^/J^cr,, 



-when ^ 0, i.e., at x = x^, and x = xr, this 



term needs to be symmetrized to be hermitian. 

Finally, the Zeeman term in the normal region is just 



Hz 



r,(T 



(A8) 



The Zeeman interaction at the leads y - {L,R} have the same 
form but in terms of the operators d^ ^ and d\.jj with summa- 
tions taken for r e Ry and replaced by gs . 



a. Conductance and Polarization 

For the normal case (Aq = 0) we compute the transport prop- 
erties in the linear response regime. The conductance follows 
from 



G = yTr[r^G''(£F)r*G"(£F)] , 



(A9) 



where C^^Xe) is the retarded (advanced) matrix propagator, 
with elements GJ^^j.,^,(fi) given by the propagator from site 
r' and spin cr' to site r and spin cr, T^^^,^, =i(S^-Spro-,r'cr' 
where Sy"^ is the retarded (advanced) self-energy due to the 
y lead and Ep is the Fermi energy. 
The polarization is 



(AlO) 



with Gcr'cr the spin-resolved conductances, i.e., the contribu- 
tions due to the electrons that are injected from lead L with 
spin cr and collected at lead R with spin cr'. As the current flow 
is in the x-direction, the relevant polarization occurs for the 
y-axis — ^recall that the Rashba SOC is proportional to p^cr,,. 
Thus, Eq.( |A10| i refers to the spin-resolved conductances with 
cr = {f , 1} along the y-axis. The polarizations along the x and 
z spin-axis are zero. 



b. Josephson current 

The Josephson current ( [T2] i, being proportional to ^ ||^, fol- 
lows from i/ftf , as it is the only term in H that depends on (p. 



2e 

/ = — i 

n 



z 

\cr,reXg 



— —e'-(c\d^^„f,} —e^-(d' ~,c^t-} + 



A* 



^^e'i0 c ,) 



(All) 
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This quantity is computed using normal propagators between 
the sites in the Xr layer and sites in the first superconduct- 
ing layer of the R lead. As Aq ^ the calculation of the re- 
quired propagators requires to work in the Nambu space since 
the equations of motion for normal and anomalous propaga- 
tors are coupled. We write all the coupled equations using the 
standard method of Ref.|44]and obtain the normal propagators 
G_„, ,(e). 



where the fermionic operators f are type c, or d according to 
the region in which they operate. 



c. Density of states, spin and current densities 

After solving the equations of motion for the propagators, 
all the densities are obtained from,"*"* 



Pra,r'Ae).Me)ds, (A12a) 

CO 

p^^yAe) = ^ {G'r^,r'A>s) " Gl^y^s)) , (A12b) 



When only the integrated quantity is important, for instance 
when studying the total Josephson current as a function of an 
external parameter, the integral in Eq.( |A12a] i can be computed 
in the complex plane by using the Residue Theorem. In this 
case, the number of integration points required for conver- 
gence is reduced by a factor of 10 as we avoid the problem 
of integrating the chain of Dirac-^s found in the real axis for 
|e| < Aq. 
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